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We propose an effective Hamiltonian approach to investigate decoherence of a quantum system in a 
non-Markovian reservoir, naturally imposing the complete positivity on the reduced dynamics of the 
system. The formalism is based on the notion of an effective reservoir, i.e., certain collective degrees 
of freedom in the reservoir that are responsible for the decoherence. As examples for completely 
positive decoherence, we present three typical decoherence processes for a qubit such as dephasing, 
depolarizing, and amplitude-damping. The effects of the non-Markovian decoherence are compared 
to the Markovian decoherence. 
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I. INTRODUCTION 

An open quantum system has been extensively studied 
for potential applications to quantum computation and 
quantum information processing 0, Qj S 3- Through 
the interaction with a reservoir a quantum system loses 
its coherent information. This is a so-called decoherence 
process 111. Decoherence has been regarded as a critical 
obstacle against quantum information processing. How- 
ever, as the notion of "reservoir engineering" has been 
suggested for laser cooled trapped ions @| , the studies on 
the decoherence have moved to control it in laboratories 
instead of suffering from it 0. Many of these studies 
have been based on Markovian reservoirs 1^. A Marko- 
vian reservoir is characterized by two essential properties: 
a) weak coupling with a system by which Born approx- 
imation is validated and b) rapid relaxation time such 
that the information of the system is diffused over the 
reservoir in a rather short time compared to the time 
scale in which the system changes. In the perspective of 
engineering a reservoir, it is desirable to study the deco- 
herence caused through a non-Markovian reservoir. It is 
strongly believed that a solid-state realization of quan- 
tum information processing may be extremely useful . 
The decoherence phenomena for solid state systems in- 
cluding a photonic band gap material and a quantum 
dot cannot always be understood under the Markovian 
assumption [IQj . The evolution of a single- mode atomic 
cavity in a non-Markovian reservoir has been studied for 
an atom laser pTf. 

A non-Markovian reservoir has "memory effects" such 
that it preserves the coherent information of the system 
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within its relaxation time. In order to analyze how many 
operations can be performed by preserving the coher- 
ent information of the system, one may consider three 
types of characteristic times: operation time Tq for a sin- 
gle quantum operation on the system, decoherence time 
T(i due to the interaction between the system and the 
reservoir, and the relaxation time within the reser- 
voir. When a quantum system decoheres in a Marko- 
vian reservoir with <^ Td, the time correlation of 
fluctuations in the reservoir is neglected. However, in 
a non-Markovian reservoir with comparable to r^, it is 
expected that the time correlation between fluctuations 
during the time interval At < Tr raises a correlated in- 
fluence on the system. Once perturbed by the system, 
the reservoir may memorize a part of the system's infor- 
mation during Tr . This memorized information 
will be fed back to the system at another perturbation 
within Tj.. It is expected that the correlated influence 
significantly suppresses the decoherence of the system 
|l4| . An extremal example for the memory effects by 
a non-Markovian reservoir is that two atoms, which are 
prepared initially in a product states, come to be in an 
entangled state through the interaction with a common 
thermal reservoir [Uli^. InRef. [13, the reservoir con- 
sists of a single-mode field with no relaxation, Tr — > oo. 
It is desirable to study decoherence processes by bridging 
the gap between the two limiting cases of Tr/r^ — > and 

Tr/Td OO. 

Dynamics of an open system may be obtained from a 
quantum Liouville equation of the total system that con- 
sists of the system and the reservoir 5, 8]. The reduced 
density operator of the system needs to be positive at all 
times. The positivity does not need to be altered by the 
presence of any other system, in other words, the reduced 
density operator of the system needs to be completely 
positive if the system and the reservoir are initially un- 
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correlated 0, 0, 0| . The complete positivity was in- 
tensively discussed for a Markovian reservoir j23| . For a 
non-Markovian reservoir, it was discussed for a master 
equation derived by using an operator sum representa- 
tion with coarse-graining approximation 2 Ij . A random 
telegraph signal model was studied for a system in a non- 
Markovian reservoir However the complete positiv- 
ity of the system density operator in a non-Markovian 
reservoir is still an important problem to be unraveled. 

In this paper, we propose an effective Hamiltonian ap- 
proach to investigate the decoherence of a quantum sys- 
tem in a non-Markovian reservoir. Instead of deriving a 
non-Markovian master equation, the formalism is based 
on the notion of an effective reservoir, i.e., certain collec- 
tive degrees of freedom in the reservoir that is responsible 
for the decoherence. We show that the reduced dynamics 
of the open system is completely positive and the complete 
positivity is naturally imposed from the effective Hamil- 
tonian approach. As examples for completely positive 
decoherence, we present three typical decoherence pro- 
cesses for a qubit such as dephasing, depolarizing, and 
amplitude-damping. The effects of the non-Markovian 
decoherence are compared to the Markovian decoherence. 



II. EFFECTIVE ENVIRONMENT 

Dynamics of an open quantum system has commonly 
been investigated by a master equation for the den- 
sity operator or a Fokker-Planck equation for the quasi- 
probability function such as a Wigner function 8]. It is 
however difficult to derive a completely positive master 
equation for a system interacting with a general non- 
Markovian reservoir. We suggest an effective Hamilto- 
nian approach by introducing the notion of "effective 
variables" in the environment. It will be shown that 
the large (normally infinite) number of environmental 
degrees of freedom, which we call "environmental vari- 
ables" , is reductive into a small number of effective en- 
vironmental variables in the sense that both cases result 
in the same master (or Fokker-Planck) equation for the 
system. For instance, a Markovian thermal environment 
is reductive into a collective single-mode boson field in a 
thermal state p^ . 

The property of being reductive into the small effec- 
tive environmental variables, that leads to the correct 
and equivalent description of the reduced dynamics for 
a system, is predicted by two observations. First, at a 
given time r, the density operator p of the system S can 
always be purified to a pure state \ip) of a larger compos- 
ite system consisting of the system S and an ancillary 
system R, such that p is obtained by tracing the pure 
state over the ancillary system, i.e., p — Trfl|V')si?(V'l- 
For the purification, an ancillary system is required to 
have its Hilbert space larger than or equal to that of the 
system [l|: For instance, a qubit may suffice for the an- 
cillary system if the system is a qubit. The composite 
system S + R provides all physical descriptions relating 



to the system and the ancillary system is called effective 
variable(s), which is in effect certain collective degrees of 
freedom, of the environment at the given time r. Sec- 
ondly, the effective variables may be dynamic over the 
whole environmental variables as the interaction time t 
passes. On the other hand, the dynamics of the effective 
variables can be absorbed by time-dependent coupling 
constants while keeping the effective variables station- 
ary. The set of the time-dependent coupling constants 
and the effective variables contains all the information 
of the environment that governs the decoherence of the 
system. We call this set an effective environment. 

Suppose Hint (t) is an interaction Hamiltonian between 
a system S and an effective environment R. We write 
the interaction Hamiltonian Hint{T), in the interaction 
picture, in the form of 

^mtW -I]^a®^a(T) (1) 

a 

where Sa is a Hermitian operator for the system S and 
Dair) — J2i3 ■^ai3{T)Rp with Hermitian operator Rp is 
for the effective environment R. The unit of fi. = 1 is 
used throughout the paper. Here, the time-dependent 
coupling constant Xaf3{T) is a Hermitian matrix. The 
composite system S + R is assumed initially in a product 
state, 

Pt(0) = /55(0)®/5fl(0). (2) 

The effective environment R consists of collective degrees 
of freedom in the environment E Its initial state 

Pr{0) is not necessarily an thermal state even when the 
real environment is thermal. The Hilbert space of the ef- 
fective environment and the time-dependent interaction 
Hamiltonian Hint{T) are determined such that they re- 
sult in the correct reduced dynamics of the system. They 
may be determined from first principle calculation such 
as the time-convolutionless projection-operator method 
[l^ . The use of the time-convolutionless projection- 
operator method is presented in Sec. III. 

Once the effective environment is determined, one may 
consider and solve the quantum Liouville equation for the 
composite system, 

^PT(r)^-i[H,^t{T),pT{T)]. (3) 
The evolution operator is now given as 

U,nt{T) = Texp(^-Z^ dT'H,^t{T')^ (4) 

where T is the time-ordering operator. In most cases, it is 
convenient to let \ai3{T) = \{T)ga(} or Da{T) = X(T)Da 
so that Uintir) = exp[— iA(r) J2a Sa'^Da] where A(r) = 
/J" dr' X{t'). This simpler form of Uint{T) has advantages 
for later discussions from the calculational point of view. 
The total system is described by the density operator 

PT{r) = UUT)pTmUr), (5) 
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and the reduced density operator of S is given by par- 
tially tracing the total density operator over R, 

Ps{t) = TriiPTir). (6) 

Letting {\n)} be the orthonormal basis set of R that di- 
agonalizes p/?,(0) = STiPnl"')fl'("'l' reduced dynamics 
of the system is described by a Kraus representation (or 
operator sum representation) (ISi] : 

ps{r) = S{t)ps{0) = Kn,n{r)ps{0)Kl„,{T) (7) 

n,m 

where Knmir) = R{n\Uint{T)\m) R^/p^^. Note that the 
superoperator 5 is a linear operator in the Hilbert- 
Schmidt space of density operators. Calligraphic letters 
are used for superoperators on the Hilbert-Schmidt space 
throughout the paper. It is remarkable that the existence 
of the Kraus representation {-/^„m (■'')} directly implies 
the complete positivity of the evolution superoperator 
S{t) [ig. This fact was already guaranteed since the 
Hamiltonian formalism was adopted in the present ap- 
proach of effective environment. This approach is appli- 
cable to the system decohering in a non-Markovian as 
well as a Markovian environment. 

The present approach of effective environment has sev- 
eral advantages in describing the reduced dynamics of a 
system. Firstly, it provides the correct and equivalent de- 
scription, similar to the master equation. Secondly and 
most importantly, it guarantees the complete positivity 
in a general environment like a non-Markovian environ- 
ment. Thirdly, it enables the analysis of the structure 
for effective variables, in a given environment, that are 
directly responsible for the decoherence. This analysis 
will be given in Sec. IIVI Lastly, one may simulate the 
decoherence in an experiment by introducing such an ef- 
fective environment and controlling the coupling to the 
given system. 



III. DETERMINATION OF EFFECTIVE 
ENVIRONMENT 

In the previous section we showed that the the effec- 
tive environment approach provides the description of 
the reduced dynamics of a quantum system, which de- 
coheres through an environment, provided the Hilbert 
space of the effective environment, its initial state Pr{0) 
and the time-dependent interaction Hamiltonian Hint (t) 
are determined. In this section we suggest a scheme 
which determines them by a first principle theory, in 
particular, the time-convolutionless projection-operator 
method 0, HH which has been employed to study 
non-Markovian environments for quantum information 
processing [m]. 



A. Time-convolutionless projection-operator 

scheme 

Consider the total system of a quantum system and an 
environment. The total system is assumed to have the 
Hamiltonian, 

Ht — Hs + He + Hint, (8) 

where Hg, He, and Hint a-i'e the Hamiltonians for the 
system, the environment, and their interaction respec- 
tively. It is straightforward to generalize that Hs may 
contain time-dependent external fields to control a quan- 
tum operation on the system. For a typical form of the 
interaction Hamiltonian, we consider a Caldeira-Leggett- 
type model given by 

H^nt = ^ 5a ® (9) 

a 

where Sa is a Hermitian operator acting on the system 
and Ea = X]fc(5ctfe^fc ~^ 9ak^k) is a fluctuating boson field 
due to perturbation arising from the system. Here, 
{e\.) is an annihilation (creation) operator for a boson 
mode k. The unperturbed boson fields of the environ- 
ment are governed by He = X^/t '^ke\ek- The set of the 
operators {Ea} describes various decoherence processes. 

The quantum Liouville equation for the density oper- 
ator pt of the total system is given by 

■^Pt{t) = ~i[HT,PT{T)] = -iCrPTir), (10) 

where Ct = Cs + Ce + ^int is the Liouville operator. 
Here, the symbol tilde is given to indicate a density op- 
erator in the Schodinger picture. Each Liouville opera- 
tor £ is a superoperator on the Hilbert-Schmidt space 
of density operators. The Liouville operators are in one- 
to-one correspondence to the Hamiltonians of the same 
subscriptions. Before the system starts to decohere, the 
environment is assumed initially in a thermal state 

= ^ exp (-iJ^/fesT) (11) 

where Z = Trexp(— iJ^/fc^T), T is the temperature 
and fcs the Boltzmann constant. The assumption may 
be released to all the time-independent states that com- 
mute with the non-interacting Hamiltonian of the envi- 
ronment, i.e., CePe = 0. 

In order to derive and solve a reduced equation for the 
system alone, we employ a projection-operator method 
[27|. Time-independent projection operators V and Q 
are defined as 

VX^PeTte{X), Q = 1-V, (12) 

for a dynamical variable X on the total system. Here Tr e 
indicates a partial trace over the environment. The re- 
duced density operator of the system is given by ps (t) = 
TtePt{t) = TTEiVpTir)]. 
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The quantum Liouville equation (|10|l can be decom- 
posed into two coupled equations for Vpr and QpT re- 
spectively by applying the projection operators. For the 
total system decoupled at r = 0, the solution for Qpx is 
substituted for the equation of VpT with an ansatz 



(13) 



The ansatz is introduced to ignore the renormalization 
of the unperturbed energy |23| (otherwise it raises the 
Lamb shift '28]). Represented in the interaction picture, 
the time-convolutionless master equation for the system 
is given ^3 by 



dr 



Ps{t) ^C{t)ps{t), 



where Ps{t) = exp{iTCs)ps{''') is the reduced density 
operator in the interaction picture and C(t) is the gener- 
alized collision operator. 

Let us consider a weak-coupling approximation (or a 



short-time limit) up to the second-order (Hi- 
collision operator C(t) can now be written as 



The 



Cir)psir) = J2 f dT'{xo.0{T-T')[Sp{T')ps{T),S„{T)] 

a/3 ^ 

+ Xo.0{t' 'T)[Sp{T),ps{r)S^{T')\} (15) 

where ^^(t) = eyiY>{iTCs)Sa and 

X«/3(-r) = Trs Ea{T)EpPE = X^a(-'r) (16) 

with £'q(t) = exp(iTCE)Ea. The function XapiT) char- 
acterizes all properties of the environment for the de- 
coherence of the system. In particular, it describes the 
time correlation between the quantum fluctuations Ea (t) 
and Ep{0) perturbed by the system at the respective 
times. The time convolutionless form of Eq. H14|l is one 
of the crucial advantages of using the projection-operator 
scheme while a time-convolution equation would be de- 
rived by a simple perturbation theory p9l | . It was shown 
[T^ that the time-convolutionless equation 1)14(1 becomes 
the Lindblad master equation in the Markov approxima- 
tion. 



B. Connection to effective environment approach 

In order to make a connection to the effective environ- 
ment scheme presented in Sec^l Eq. ((15|l will be further 
analyzed. In Eq. ((15(1 the system operators depend on the 
evolution time. The time-dependent operators can be ex- 
panded by a complete set of Hermitian operators, C, in 
the Hilbert-Schmidt space where the set C is chosen to 
include the operators Sa of the system: 



Sair) = ^ Cai3{T)Sf3. 

13=1 



(17) 



where d is the dimension of the systems Hilbert space. 
The coefficient Cq^(t) can be obtained by solving the 
Heisenberg equation, 

^S^{t)^i[Hs,L{t)]. (18) 

Eqs. ((17(1 and 1(181) give the following equation, which is 
equivalent to Eq. I(15() . 



C{r)ps{r) 



X! {lal3{T)[SpPs{T),Sa\ 
(1/3=1 

+7^a(T)['S/3,/5s(T-)5a]| 



(14) where 



dr'xapir, t') 

dT'Xal3{T',T) 



Xq/3(t,t') = ^Cya{T)x^s{T - t')csi3{t') 
7(5 



(19) 

(20) 

(21) 
(22) 



In the form of Eq. p9(l . the matrix 7^/3 (r) determines 
the master equation ((14(1 and the reduced dynamics of the 
system. The matrix "/af3 (t) may be represented by a sum 
of Hermitian and anti-Hermitian matrices, i.e., 7q/3(t) = 
lapi"^) + *7a/3 (''')• 'l'^^ Hermitian term 7^(t) involves 
in the decoherence while the anti-Hermitian term 17^^ (t) 
contributes to the Hamiltonian dynamics of the system. 
We call the Hermitian matrix 7^ (r) a decoherence rate 
matrix. Now, the system operators in Eq. 1(19(1 are time- 
independent which enables the connection to the effective 
environment approach. 

The evolution superoperator S'{t), i.e., the solution to 
the master equation l(14() . is not necessarily completely 
positive in the long-time limit. On the other hand, the 
method of the effective environment guarantees the com- 
plete positivity of the evolution superoperator 5(r) at 
all times. In the Appendix we show how to examine the 
complete positivity for a given superoperator. In the ef- 
fective environment approach one needs to determine the 
parameters of the effective environment. To determine 
them we require that in the short-time limit the mas- 
ter equation derived by the effective environment be ap- 
proximately equal to the master equation ((14(1 in conjuc- 
tion with Eq. ((19(1 . The effective environment approach 
may have different higher-order terms from the time- 
convolutionless projection-operator method and this dif- 
ference leads to the complete positivity of the reduced 
dynamics in the effective environment approach. 

Using the interaction Hamiltonian ^ and following 
the similar procedure leading to Eq. 1(141) . we can de- 
rive the master equation in the effective environment 
approach. This has to be approximate to Eq. ((14() in 
conjuction with Eq. ((19() in the short-time limit. This 
comparison leads to the following two conditions 



Tr, 



^a(T)/5fl(0) =0, 



(23) 
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which is comparable to Eq. H13|) . and 

dr'TvR \D^iT)Df3iT')pR{0)\ « 7„0(t) (24) 



where Daij) is given in Eq. The initial density 

operator Pr{Q) and the time-dependent coupling matrix 
Aa/3(''") can be chosen by varying the operators Ra so 
as to satisfy the coupled linear equations (|23|l and (|24|l . 
The Hilbert space of the effective environment needs to 
be enlarged with respect to the dimensionality unless any 
solution is found to Eqs. (|23|l and H24|l . As a result, the 
chosen set of operators {Dair)} are responsible for the 
type of decoherence which the system undergoes. In this 
sense the operator Da (t) is called a decoherence channel 
to the environment. 

In most cases, an effective environment has the same 
Hilbert-space dimensionality as the quantum system. For 
instance, this is the case for Markovian environments. 
Nonetheless, an effective environment can be of larger 
dimensionality in order to properly describe the reduced 
dynamics of a system. An effective environment may be 
classified into a "Q-qubit environment" if the effective 
environment consists of Q qubits, "M-mode bosonic en- 
vironment" if it consists of M bosonic modes, and so on. 

We shall consider the asymptotic form of the time- 
correlation function x(r) = ^7(t). The time-correlation 
function is in general complex. However, being interested 
in the decoherence, we concentrate on the real part of the 
time-correlation function. The time-correlation function 
becomes real- valued when the system Hamiltonian is re- 
defined so that it includes Lamb shifts Noting that 
the time-correlation disappears at a long time, x(''") ap- 
proaches zero in the limit of r ^ oo. Thus, x(t) is 
assumed to have the form of 



x(t) = /(r)e-^'M 



(25) 



where /(r) and gir) are polynomial (or trigonometric) 
functions and further gir) ^ oo as r — > ±c». A simple 
form of xi^) will be given by 



X(r) 



4t. 



(26) 



where k — l/r^ is a decoherence rate with the deco- 
herence time and is a relaxation time (or memory 
time) in which the injected information is diffused over 
the environment. In the Markovian limit of <ti t^, 
the time correlation function becomes a delta function 
so that x{t) ~^ k6{t)/2 and j{t) k/4. The time- 
correlation function in Eq. (|26(l may be regarded as the 
first extension from Markovian to non-Markovian deco- 
herence [T^. 



IV. QUBIT SYSTEM 

We shall present typical decoherence of a qubit system 
in Markovian and non-Markovian environments in the 



approach of an effective environment. In addition, we 
investigate the structure of an effective environment: a) 
the environmental variables that play a role in the respec- 
tive decoherence, b) the coupling constants between the 
system and the effective environment, and c) the initial 
quantum states of the effective environment. 

We shall first of all consider a dephasing process. As 
in Sec. IIIII we compare both master equations in the 
projection-operator method and the effective environ- 
ment approach for a dephasing process. It is found that 
the effective environment has one qubit and the interac- 
tion Hamiltonian is given by 



Hintir) = A(t)o-^ (g) o-^, 

27(T)exp[-4r(T)] 



A(r) 



v/1 - exp[-8r(r)] 



(27) 
(28) 



where CTq, is a Pauli spin operator and T{t) = dr'^lr'). 
The initial states of the system and the one-qubit effec- 
tive environment can be written as 

PsiO) - i(i + s.^), (29) 



MO) 



(30) 



where 1 is an identity operator, s = {sxTSy, Sz) a Bloch 
vector of the system qubit, r = {rx,ry,rz) a Bloch vec- 
tor of the one-qubit environment, and & — ((ij., CTj,, ct^). 
In connection with the time correlation function of the 
environment, the coupled linear equations (|23|l and H24|) 
are satisfied for r = {rj;,ry,0) and A(r) in Eq. in 
the short-time limit of |7(T)|r ^ I. The initial state 
of the effective environment is unpolarized along the z 
axis while it may be polarized on the x-y plane for the 
dephasing process. 

Let us consider the reduced dynamics of the system. 
Applying Eqs. and ij^HIl to Eqs. lEI)-®, we obtain 
the reduced dynamics in terms of its Bloch vector. 



sir) = (5,exp[-4r(r)],Syexph4r(T)], 



(31) 



The dynamics of the dephasing process is characterized 
by the decay function exp[— 4r(r)]. The decay function 
may be regarded as a fringe visibility on the quantum 
interference for a given time r. In order to compare the 
effect of the non-Markovian decoherence to the Marko- 
vian one, we employ the simple time-correlation function 
in Eq. (|26|l . The decay function then becomes 



exp[— 4r(r)] — exp |^'* t — Tr 



1 - e-^/^ 



} ■ (32) 



In the limit of no memory, — > 0, the dephasing process 
is Markovian with the decay function, exp[— 4r(T)] = 
exp(— Atr). When the environment keeps the information 
of the system within the memory time Tr, the loss of the 
coherent information is suppressed in comparison to the 
Markovian case as present in Fig. ^ This fact can be 
clearly seen through the decoherence rate, 



7(t) 



dr 



r(r) 



i.(f-e-^/^^) 



(33) 
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for a given memory time Tr- The decoherence rate is k/4 
if — > and it is reduced by increasing the memory 
time r,.. 




KT 



FIG. 1: Decay function, exp[— 4r(r)], with respect to the 
evolving time r. The decay function may be regarded as a 
fringe visibility for quantum interference as the time passes. 
It depends on the memory time Tr for a given decoherence 
rate n = l/rd'. a) kt^ — 0.1 (solid line), b) ktt = 1 (dotted 
line), and c) kt^ — 10 (dashed line). The case a) is already 
close to the Markovian decoherence {ktv — > 0) 

Similar analysis can be applied to depolarizing or 
amplitude-damping processes for a qubit. For this pur- 
pose, the interaction Hamiltonian is assumed to be 

gapCTa^Crp (34) 

a.j3—x.y,z 

where A(r) is given by Eq. (^5)) and gap is a real coupling 
matrix which determines the type of decoherence. The 
initial states of the system and the one-qubit effective 
environment are given by their Bloch vectors s and r, 
similar to Eqs. H29I) and H3UI) . The coupling matrix g^p 
and the initial state r of the effective environment are 
summarized in Table. ^ It is interesting to look at the 
initial state of the effective environment for each process. 
The initial state is completely random with r = (0, 0, 0) 
for the depolarizing process and it is polarized along z 
axis for the amplitude-damping process. 

The reduce dynamics of the qubit is described in terms 
of the time dependence of its Bloch vector s(r). For the 
depolarizing process, the time-dependent Bloch vector is 
given by 

Spo!(T) = sexp[-8r(r)] (35) 

For the amplitude damping process, it is given by 

[samp(T)]^ = Sj; exp[-4r(r)] 
Kmp{T)]y = Sj,exp[-4r(r)] 
Kmp{T)], = + (s, -r,)exp[-8r(T)]. (36) 

It is well known that a system slowly decoheres in a 
non-Markovian environment due to the memory effects. 



TABLE I: The coupling matrix ga/i and the initial state r 
of the effective environment for dephasing, depolarizing, and 
amplitude damping processes. 





gap 


r 


dephasing 






depolarizing 




(0,0,0) 


amplitude damping 


SalsiS/Sx + Spy) 


(0,0, r,) 



However, we stress that the result should be derived from 
the completely positive reduced dynamics and our anal- 
ysis of the effective Hamiltonian approach strongly sup- 
ports the result. 

V. REMARKS 

The effective Hamiltonian approach was proposed to 
investigate the non-Markovian decoherence of an open 
system and further to understand the characteristics of 
its environment. The formalism is based on the notion of 
an effective environment, i.e., certain collective degrees of 
freedom in the environment that are responsible for the 
decoherence. The present approach naturally imposes 
the complete positivity on the reduced dynamics for the 
system. We applied the approach to the dephasing, depo- 
larizing, and amplitude-damping processes. It was found 
that the non-Markovian environment suppresses the de- 
coherence of the qubit, due to the memory effect of the 
non-Markovian environment. 
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APPENDIX A: COMPLETE POSITIVITY OF A 
SUPEROPERATOR 

Consider a set of bounded operators acting on the state 
vector in the Hilbert space. The set of bounded operators 
forms a vector space V with the inner product defined 
by Hilbert-Schmidt norm [s^- The inner product of two 
bounded operators v and w in y is given by 

{{w\v)) EE Tt {w^ v) . (Al) 

In the analogy of the ket and bra states, \v)) is called 
a "ket" vector of an operator ii and {{v\ a "bra" vector. 
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The set of the bra vectors forms the dual vector space. 
Let {cab = be an orthonormal basis for the vector 

space with the completeness relation, 



\eab)){{eab\ =1 



(A2) 



ab 



where X is an identity superoperator, i.e., a linear map of 
the ket vector (or ket operator) to itself. Any ket vector 
\v)) is expanded as 



(A3) 



where Vab = {{eab\v)). 

A superoperator 5 is a linear map of V onto itself 



S 



(A4) 



The matrix elements of S in the orthonormal basis {iab} 
is obtained as 

Sab.cd = ((eafc|5|ecd)) = ((eab|5(ecrf))). (A5) 
The superoperator S is called "Hermitian" when 

A superoperator S is called "positive" when {{fi\S\fj)) is 
real and positive for all fj G V. 

If S is Hermitian, it has the right eigenvector \vs)) with 
the "real" eigenvalue such that 



(A7) 



and {{vs\ is the corresponding left eigenvector with the 
same eigenvalue. Further, if and only if S is positive, the 
eigenvalues are all positive. When iS is not Hermitian, 
it may have complex eigenvalues and {{vs\ is no longer 
the corresponding left eigenvector. Instead, there exists 
a left eigenvector ((ws| such that 



({Ws\S = ((-Ws\Xs 



(A8) 



where As is a complex number. The right and left eigen- 
vectors satisfy the orthogonality 

{{w,\v,^)) ^ {{w,\v,))6ss'. (A9) 
The completeness relation is given as the following 



{{Ws\Vs)) 

and S can be represented as 



I 



{{Ws\Vs)) 



(AlO) 



(All) 



Definition: Partial transposition, denoted by "#" , on 
a superoperator S is defined as 



A superoperator S is said to be Hermitian-preserving 
when S{x^) = S{xy for any x G V. The partial trans- 
posed superoperator S'^ is Hermitian if and only if S is 
Hermitian-preserving. In addition, S"^ is closely related 
to the complete positivity of S. 

Theorem: The following three conditions for a super- 
operator S are all equivalent: 

(a) S is completely positive 

(b) S has a Kraus representation, i.e., 



S{x) = Y,K^xkl 



(A13) 



(c) 5* is positive 

Proof: It is well known ,18j that (a) and (b) are equiv- 
alent. Consider (b) ^ (c). When S has a Kraus repre- 
sentation, its matrix elements are given as 



bd* 



Sab.cd = ^ ^r^A^' 

The partial transposition S'^ becomes 

C'# C \ ^ I^ib l^'^d* 

•^ab,cd — ^ac,bd — ^ -"-^ 



(A14) 



(A15) 



It is clear that for all 77 e F 



ibcd 



where = J2ab Vlb-^if- Thus 5* is positive. 

The converse (c) ^ (b) is now considered. Since S'^ 
is positive, it is also Hermitian and all the eigenvalues 
are positive: 



(A17) 



where \vi,)) are normalized eigenvectors. Since di, is pos- 
itive, we can define the following matrix 

Kf = Vd^{{eab\v.))- (A18) 
Thus the superoperator S has a Kraus representation as 

Sab,cd = St^,, = 5](i^r)(^r)*- (A19) 

The set of Kraus operators are not unique 0| and all 
equivalent sets of Kraus operators for the given super- 
operator can be generated by "unitary remixing" of the 
canonical set with the eigenvalue vector d' extended b"v 



canonical set witn tne eigenvalue vector d extended by 
Stbcd = {{^ab\S*\ecd)) = ((eacl^lebd)) = S'ac,6rf-(A12) some arbitrary number of zeros, i.e. d' = (d, 0,...,0) [3l|. 
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